Identification of plasma miR-4505, miR-4743-5p and miR-4750-3p as novel diagnostic biomarkers for coronary artery disease in patients with type 2 diabetes mellitus: a case-control study

Background Type 2 diabetes mellitus (T2DM) and coronary artery disease (CAD) are commonly coexisting clinical entities with still growing incidence worldwide. Recently, circulating microRNAs (miRNAs) have emerged as novel molecular players in cardiometabolic diseases. This study aimed to identify a specific miRNA signature as a candidate biomarker for CAD in T2DM and to delineate potential miRNA-dependent mechanisms contributing to diabetic atherosclerosis. Methods A total of 38 plasma samples from T2DM patients with and without CAD, CAD patients and healthy controls were collected for expression profiling of 2,578 miRNAs using microarrays. To investigate the regulatory role of differentially expressed (DE)-miRNA target genes, functional annotation and pathway enrichment analyses were performed utilizing multiple bioinformatics tools. Then, protein-protein interaction networks were established leveraging the STRING database in Cytoscape software, followed by cluster analysis and hub gene identification. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) was carried out for microarray data validation in the larger replication cohort of 94 participants. Receiver operating characteristic analysis was applied to evaluate the diagnostic values of miRNAs. Multivariate logistic regression analysis was used to develop miRNA-based diagnostic models. Results In the discovery stage, overexpression of hsa-miR-4505, hsa-miR-4743-5p, hsa-miR-6846-5p, and down-regulation of hsa-miR-3613-3p, hsa-miR-4668-5p, hsa-miR-4706, hsa-miR-6511b-5p, hsa-miR-6750-5p, hsa-miR-4750-3p, hsa-miR-320e, hsa-miR-4717-3p, hsa-miR-7850-5p were detected in T2DM-CAD patients. The DE-miRNA target genes were significantly enriched in calcium ion binding, regulation of actin cytoskeleton, and gene expression. hsa-miR-4505, hsa-miR-4743-5p, and hsa-miR-4750-3p were found to be involved in fatty acid metabolism, leukocyte transendothelial migration, and neurotrophin signaling pathway. Dysregulation of hsa-miR-4505, hsa-miR-4743-5p, and hsa-miR-4750-3p in T2DM-CAD patients compared with T2DM subjects and controls (all p < 0.001) was further confirmed by RT-qPCR. All validated miRNAs demonstrated good discriminatory values for T2DM-CAD (AUC = 0.833–0.876). The best performance in detecting CAD in T2DM was achieved for a combination of three miRNAs (AUC = 0.959, 100% sensitivity, 86.67% specificity). Conclusions Our study revealed a unique profile of plasma-derived miRNAs in T2DM patients with CAD. Potential miRNA-regulated pathways were also identified, exploring the underlying pathogenesis of CAD in T2DM. We developed a specific three-miRNA panel of hsa-miR-4505, hsa-miR-4743-5p and hsa-miR-4750-3p, that could serve as a novel non-invasive biomarker for CAD in patients with T2DM. Supplementary Information The online version contains supplementary material available at 10.1186/s12933-024-02374-0.


Background
Type 2 diabetes mellitus (T2DM) is a chronic metabolic disease characterized by insulin resistance and hyperglycemia, which play a crucial role in accelerating endothelial damage and the onset of cardiovascular complications, making individuals with diabetes two-to four-times more likely to develop coronary artery disease (CAD) than those without diabetes [1].Additionally, the mortality of patients with T2DM is considerably higher in the presence of cardiovascular diseases, accounting for more than 50% of all deaths in this population [2].The phenomenon of CAD in T2DM is multifactorial and involves an interplay of genetic, epigenetic and environmental factors, but only 25% of them have been identified so far [3,4].It is well known that T2DM-associated CAD has an earlier onset, faster progression, and greater extent compared to non-diabetic subjects [5].What is more, approximately one-third of asymptomatic individuals with T2DM were found to have significant CAD with a subsequent increased risk of major adverse cardiac events [6].Hence, early detection of CAD in patients with T2DM is becoming an important public health concern.
Currently, invasive coronary angiography remains the gold standard in the diagnosis of CAD, but it is reserved for patients whose clinical risk is assessed as high or when stress testing indicates significant ischemic burden [7,8].Moreover, the routine use of this technique is limited due to its invasiveness and potential side effects of ionizing radiation, contrast agent, which may generate various procedure-related and patient-dependent complications [9,10].Therefore, it is highly necessary to find alternative non-invasive tools that can provide additive predictive information to improve risk stratification of CAD in patients with T2DM, thus avoiding unnecessary invasive testing and actively preventing the development of diabetic atherosclerosis.
MicroRNAs (miRNAs, miRs) are a family of small (18-25 nucleotides) non-coding RNA molecules that modulate gene expression at the post-transcriptional level by inhibiting the translation of target mRNAs or inducing mRNA degradation [11].Each specific miRNA has the ability to target multiple transcripts, while a single mRNA can be regulated by a set of miRNAs, allowing precise control of a wide range of biological processes [11].That is why, the establishment of a specific miRNA panel may better reflect the dynamic changes of these molecules in response to disease states than a single miRNA [12].As the number of miRNAs continues to increase, simultaneous expression profiling of hundreds of miRNAs in a single experiment appears to be indispensable to evaluate novel miRNA-mRNA interactions and clinical applications of miRNAs, particularly those not previously screened for disease [12].Besides, it enables to avoid miRNA pre-selection bias based solely on literature overview [12].In contrast to other RNAs, miRNAs have proven to be extremely stable in circulation and resistant to enzymatic degradation by RNase and deleterious external conditions [13].miRNAs are also readily detectable in a variety of body fluids, including serum and plasma, highlighting their potential role as disease biomarkers [13].
Since bi-directional links between T2DM and CAD exist, the discovery of miRNAs as essential epigenetic regulators of common pathological events in cardiovascular and metabolic diseases have attracted a tremendous interest [14].Altered expression of endothelial-specific miR-92a, miR-126, miR-210 along with miR-342 and miR-450 is known to promote hyperglycemia-induced oxidative stress and inflammation accelerating endothelial dysfunction [15][16][17][18].Moreover, miRNAs have been reported to trigger critical processes involved in diabetic atherosclerosis, including monocyte activation, leukocyte transendothelial migration, vascular smooth muscle cell (VSMC) proliferation/migration, and platelet reactivity [14,19].Circulating miRNAs may strictly correspond to tissue injury, which raises endothelial-enriched miR-17, miR-92a, miR-126, cardiac-specific miR-1, miR-133, miR-499, and VSMC-enriched miR-145 as fingerprints for diabetic atherosclerosis [14].Although, the diverse role of miRNAs in the etiology of T2DM or CAD has already been extensively explored, studies evaluating the utility of miRNA-based biomarkers in coexisting CAD and T2DM are scarce [14].
The present study aimed to investigate the differential expression of circulating miRNAs in plasma samples from T2DM patients with and without CAD, and to identify a specific miRNA signature that could serve as a potential biomarker for CAD in T2DM.In addition, target annotation and functional enrichment analyses were performed to delineate potential miRNA-dependent mechanisms contributing to the development of CAD in T2DM individuals.

Study design and patients
This study was designed as an observational, case-control study using a two-stage approach, including primary screening and ensuing validation stage.All recruited participants were Caucasian individuals hospitalized in the Department of Endocrinology, Diabetology and Metabolic Diseases at the Medical University of Lublin, Poland, between October 2020 and December 2022.In the discovery phase, 38 patients (T2DM with CAD group, n = 12; T2DM group, n = 12; CAD group, n = 8; control group, n = 6), matched for age, sex, and body mass index (BMI), were selected for the miRNA expression profiling with microarrays.Candidate miRNAs were subsequently validated in the independent replication cohort of 94 age-, sex-, and BMI-matched patients (T2DM with CAD group, n = 30; T2DM group, n = 30; CAD group, n = 16; control group, n = 18) using reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR).The overall design and key stages of the study are outlined in the flow chart (Fig. 1).
The diagnosis of T2DM was established according to the American Diabetes Association criteria [20].Patients with T2DM and/or CAD underwent invasive coronary Fig. 1 Flow chart of the study.T2DM-CAD, type 2 diabetes mellitus with coronary artery disease; T2DM, type 2 diabetes mellitus; CAD, coronary artery disease; DE-miRNA, differentially expressed miRNA; FC, fold change; FDR, false discovery rate; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; DIANA-miRPath, DNA Intelligent Analysis-miRPath; PPI, protein-protein interaction; STRING, Search Tool for the Retrieval of Interacting Genes; MCODE, Molecular Complex Detection; RT-qPCR, reverse transcription quantitative real-time polymerase chain reaction; ROC, receiver operating characteristic angiography which was assessed by two independent interventional cardiologists, and CAD was defined as at least one major epicardial vessel with > 50% stenosis, while those without any stenosis or luminal irregularities in any of the epicardial arteries were classified as CAD negative ones [21].The control subjects were randomly selected from healthy volunteers who had normal glucose tolerance and no evidence of CAD.
Detailed information about demographic data, medical history, current medication, alcohol consumption, and smoking habits were obtained from each participant.The study was limited to a population aged 45-65 years.All individuals in the study groups were on a stable hypoglycemic, antihypertensive, lipid-lowering, and/ or antiplatelet treatment regimen for at least 3 months prior to participation in the study.The exclusion criteria for the study encompassed the following: (1) other types of diabetes mellitus, (2) micro-and other macrovascular diabetic complications (retinopathy, neuropathy, nephropathy, peripheral artery disease, stroke, cerebrovascular disease) or acute diabetic complications (diabetic ketoacidosis and coma, hyperglycemic hyperosmolar state, lactic acidosis), (3) autoimmune diseases, (4) acute and chronic inflammatory diseases, (5) human immunodeficiency virus or hepatitis C virus infection, (6) liver dysfunction (aspartate aminotransferase or alanine aminotransferase levels exceeding 2 times the upper limit of the normal range), (7) renal impairment with estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73m 2 , (8) cancers, (9) previous myocardial infarction, percutaneous coronary intervention and/or coronary artery bypass grafting, (10) heart failure, (11) cardiac arrhythmias, (12) moderate to severe valvular heart disease, congenital heart disease, cardiomyopathy, (13) surgery or trauma in the past 3 months, (14) the use of drugs with the proven effect on blood counts, including glucocorticoids, in the past 6 months, (15) a daily alcohol consumption greater than or equal to 30 g for men and 20 g for women, (16) smoking.
The study protocol was approved by the Bioethics Committee of the Medical University of Lublin, Poland (No. KE-0254/198/2020) and the research was conducted in accordance with Good Clinical Practice (Declaration of Helsinki of 1975, revised in 2013).Written informed consent was obtained from all participants prior to enrollment in the study.

Anthropometric measurements
Anthropometric parameters such as body weight and height, waist and hip circumferences were measured (to the nearest 0.1 kg and 0.1 cm) using standardized equipment and methodology.BMI was calculated as weight in kilograms divided by height in meters squared (expressed in units of kg/m 2 ), and waist-to-hip ratio (WHR) as the ratio of waist circumference to hip circumference.Systolic (SBP) and diastolic (DBP) blood pressure were taken in three sets using an oscillometric method in a sitting position after a minimum of 5 min of quiet rest, and calculated as the average of the second and third measurements.A twelve lead electrocardiography was performed on all subjects and was carefully reviewed by a cardiologist for abnormalities.

Sample collection
Whole blood (10 mL) was collected by venipuncture after a 12-hour overnight fast (between 7:30 and 8:30 a.m.) into plastic tubes (S-Monovette, Sarstedt, Nümbrecht, Germany) and processed within 30 min to obtain plasma and serum for routine laboratory analyses.Biochemical measurements, including plasma glucose, glycated hemoglobin A1c (HbA1c) and fibrinogen, serum triglycerides (TG), total cholesterol (TC), high-density lipoprotein (HDL)-C, uric acid, creatinine, high-sensitivity C-reactive protein (hs-CRP), and total homocysteine concentrations were performed by automated standard methods in a centralized laboratory using Atellica CH 930 analyzer (Siemens Healthineers, Erlangen, Germany).The concentration of low-density lipoprotein (LDL)-C was calculated according to Friedewald equation.eGFR was calculated according to Modification of Diet in Renal Disease (MDRD) formula.Complete blood counts were analyzed using an automated hematology analyzer (Yumizen H500 Analyzer, Horiba Medical, Northampton, UK).
Five millilitres of whole blood for miRNA analysis were collected in ethylenediaminetetraacetic acid (EDTA)containing plastic tubes (S-Monovette, Sarstedt, Nümbrecht, Germany) at the same time and under the same conditions as the samples for biochemical testing.The samples were stored on ice and processed immediately after the draw using two-step centrifugation protocol: centrifugation at 3,000× g for 10 min at 4 °C to remove blood cells, and then centrifugation of the isolated supernatant at 16,000× g for 10 min at 4 °C to remove additional cellular nucleic acids attached to cell debris and to obtain platelet-poor plasma.Subsequently, the plasma was transferred to RNase/DNase-free tubes (DNA LoBind Tubes, Eppendorf, Hamburg, Germany) and stored at − 80 °C, in the absence of freeze-thaw cycles, until further analysis.

Total RNA isolation
The total RNA, including miRNA, was extracted from 200 µl of plasma using the miRNeasy Serum/Plasma Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions.The RNA concentration was determined using a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).For all analyzed samples, the A260/A280 purity ratio was between 1.8 and 2.0.In addition, the quality and integrity of RNA was assessed with an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA) and the RNA 6000 Pico Kit according to the manufacturer's protocol.The isolated RNA was stored at − 80 °C for further analysis.

miRNA expression profiling using microarrays
The miRNA expression profile was determined using the microarray system with the GeneChip™ miRNA 4.0 Array chip (Affymetrix, Santa Clara, CA, USA).Briefly, the total RNA (400 ng) was labeled using the FlashTag™ Biotin HSR RNA Labeling Kit (Affymetrix, Santa Clara, CA, USA) following the manufacturer's recommendations.The labeled RNA was hybridized to gene chips in a GeneChip™ Hybridization Oven 645 (Affymetrix, Santa Clara, CA, USA) at 48 °C for 16 h.The gene chips were washed and stained using the GeneChip Hybridization Wash and Stain Kit (Affymetrix, Santa Clara, CA, USA), and then scanned with a GeneChip™ Scanner 3000 7G (Affymetrix, Santa Clara, CA, USA).
The raw intensities in a CEL format were imported into the Transcriptome Analysis Console (TAC) software version 4.0.1 (Thermo Fisher Scientific, Waltham, MA, USA), normalized using the Robust MultiArray Average (RMA) method, and processed by the Detection Above Background (DABG) algorithm to determine whether a given probeset is reasonably expressed (p-value < 0.05).Then, the data were log2-transformed for further analysis to determine differentially expressed miRNAs (DE-miR-NAs) between patients with T2DM, T2DM-CAD, CAD and controls.miRNAs that showed fold change (FC) of >|1.5| between the study and control groups, with a p-value < 0.05 and an adjusted p-value by the Benjamini-Hochberg correction for multiple hypothesis testing (False Discovery Rate, FDR value) ≤ 0.05, were considered statistically significant and differentially expressed.

Validation of microarray results by reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR)
To confirm the expression of selected miRNAs identified by the microarray analysis, RT-qPCR was performed.In the first step, 1 µg of total RNA was reverse transcribed into first-strand cDNA using the TaqMan™ miRNA Reverse Transcription Kit (Thermo Fisher Scientific, Waltham, MA, USA) with the addition of an RNase inhibitor (Thermo Fisher Scientific, Waltham, MA, USA).The reaction mixture was then incubated at 16 °C for 30 min, at 42 °C for 30 min and at 85 °C for 5 min in a Veriti™ Dx 96-well Thermal Cycler (Applied Biosystems, Thermo Fisher Scientific, Foster City, CA, USA), and then held at 4 °C.cDNA samples were stored at − 20 °C until analysis.
The reactions were incubated in a 96-well MicroAmp Fast Optical 0.1 mL reaction plate (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA) at 50 °C for 2 min and at 95 °C for 10 min, followed by 40 cycles of PCR amplification at 95 °C for 15 s, and elongation at 60 °C for 1 min.All the PCR reactions were run in triplicate.

Assay name
Assay ID miRBase ID 1 miRBase Accession Number 1   Primer sequence (5' to 3') According to miRBase 22 (https://www.mirbase.org/) 2 NCBI Accession Number The relative expression (RQ) of each individual miRNA was calculated using the 2 −ΔΔCt method, and logarithmically transformed (log 10 RQ) to reduce skewness [22].To normalize the expression levels of candidate miR-NAs, the cycle threshold (Ct) value was determined for each sample against the reference gene, small nuclear RNA (snRNA) U6.In the study, RNU6B (RNU6-6P) and RNU48 (SNORD48) were also considered as endogenous controls due to their relatively constant and highly abundant expression across human tissues and cell line types, and a wide use in various fields, including diabetes and coronary artery disease research, but U6 showed the most stable expression in the tested material, and therefore proved to be the best control.The data were analyzed with Expression Suite Software version 1.0.3.(Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA) with the automatic Ct setting for assigning baseline and threshold for Ct determination.

Functional enrichment analysis of differentially expressed miRNAs
The web-based DNA Intelligent Analysis (DIANA)-miRPath v3.0 tool (https://dianalab.e-ce.uth.gr/html/mirpathv3/index.php?r=mirpath, accessed on 31 May 2023) was used to perform functional enrichment analysis for each of the DE-miRNAs in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) terms [23].The target genes of DE-miR-NAs were obtained by combining available in silico predicted targets from the DIANA-microT-CDS v5.0 and/ or TargetScan v6.2 databases and high quality, manually curated and experimentally validated targets from the DIANA-TarBase v7.0 database [23].The p-values were calculated with Fisher's exact test, and adjusted using FDR correction for multiple comparisons.The cut-off criterion was a p-value limited to 0.05.Visualization of the miRNA regulatory network with appropriate KEGG and/or GO annotations was carried out using Cytoscape v3.10.0 software (https://cytoscape.org/,accessed on 15 June 2023) [24].
Subsequently, an enrichment analysis was performed using the miRNet 2.0 tool (https://www.mirnet.ca/miR-Net/home.xhtml,accessed on 25 June 2023), where a list of selected miRNAs and a list of genes with which they interact were introduced [26].KEGG and Reactome pathways and GO categories were searched to understand the biological relevance of DE-miRNA target genes in T2DM-CAD.The terms with an adjusted p-value < 0.05 was regarded as statistically significant.The plot visualizing the enrichment analysis was generated using the ggplot2 3.3.0package in the R environment.

Protein-protein interaction networks and sub-networks construction
To establish the regulatory role of validated DE-miRNA target genes in T2DM-CAD, the Search Tool for the Retrieval of Interacting Genes (STRING) App v2.0.1 (https://apps.cytoscape.org/apps/stringapp,accessed on 25 July 2023) in Cytoscape v3.10.0 software was used to construct the protein-protein interaction (PPI) networks [27].Interactions between proteins encoded by coexpressed genes, including direct (physical) and indirect (functional) associations, were obtained from different sources, and only interactions with a combined confidence score greater than 0.4 were retained.
Then, the Molecular Complex Detection (MCODE) v2.0.3 (https://apps.cytoscape.org/apps/mcode,accessed on 25 July 2023) plugin of Cytoscape v3.10.0 tool was used to gain insights into possible closely connected gene modules (clusters, sub-networks) within the PPI networks and perform an enrichment analysis [28].The networks were clustered based on inflation parameters according to the Markov Clustering (MCL) algorithm, and the 6 top-scored clusters separately for up-and down-regulated miRNAs were analyzed.The selection criteria were as follows: cluster finding = hair-cut, node score cut-off = 0.2, degree cut-off = 2, k-core filter of 2 and maximum depth = 100.In the networks, each node represents a protein produced by a single, protein-coding gene locus, and edges indicate both functional and physical associations of the proteins.

Hub gene identification
To identify the key genes (hub genes) in the PPI networks, the CytoHubba plugin (https://apps.cytoscape.org/apps/cytohubba, accessed on 25 July 2023) of Cytoscape v3.10.0 software was applied [29].The top 10 ranked hub genes were extracted using the BottleNeck algorithm according to their betweenness centrality.Next, an enrichment analysis was conducted to forecast the function of identified hub genes of the T2DM-CAD gene regulatory networks.

Statistical analysis
Sample size calculation was based on similar previous experiments and pilot data on microarray miRNA expression profiling in patients with T2DM and CAD [30][31][32][33].We used the MD Anderson sample size calculator (http://bioinformatics.mdanderson.org/Microarray-SampleSize)to estimate the minimum number of samples per experimental groups (T2DM-CAD and T2DM) to detect 1.5 fold differences in miRNA expression levels between groups at the true positive detection power (1-β) of 80% [34][35][36].The sample size of 12 patients for each of the 2 groups was calculated to provide 80% power with a significance level (α) set at 5% and a standard deviation (SD) of 0.5.
Categorical variables are expressed as numbers with percentages (%).Continuous data are presented as mean ± SD or median with interquartile range (IQR), depending on the normality of distribution assessed by the Shapiro-Wilk test.Parametric tests were used for variables with normal distributions.Student's t-test was used to compare differences between two groups, whereas one-way analysis of variance (ANOVA) test was used for more than 2 groups with Tukey's post hoc test for unequal n.Data conforming to non-normal distribution were analyzed using the Mann-Whitney U test (for comparisons of two groups) and the Kruskal-Wallis one-way analysis of variance by ranks followed by Dunn's test (for comparisons of more than two groups).Chi-Squared test was performed to determine sex distribution between groups.
To estimate the diagnostic significance of particular DE-miRNAs in detecting CAD in T2DM patients, receiver operating characteristic (ROC) analysis was performed, and the area under the curve (AUC) was calculated.Multivariate logistic regression analysis using the generalized linear model (GLM) function included in R was used to train miRNA diagnostic models.The expression value of identified miRNAs was considered as predictive variables, and the sample type (T2DM-CAD or T2DM) was considered as a binary response variable.To assess the diagnostic accuracy of each multi-miRNA panel, ROC curves and AUC were calculated by the pROC package.The cut-off points for ROC curves were calculated by the Youden index.
All p-values are two-tailed, and p-values < 0.05 were considered as statistically significant.The statistical analysis was performed with STATISTICA 13.

Baseline characteristics of study and control groups
A total of 132 patients meeting eligibility criteria were prospectively recruited for the miRNA profiling study (n = 38) and the validation cohort (n = 94).All participants in both cohorts were divided into four groups (T2DM-CAD, T2DM, CAD, controls) matched for age, sex, and BMI.As expected, patients with T2DM had significantly higher HbA1c levels and fasting plasma glucose concentrations compared to non-T2DM individuals (p < 0.001).None of the study groups differed in a series of clinical and biochemical parameters (SBP, DPB, TC, TG, LDL-C, hs-CRP, platelets, fibrinogen, uric acid, and renal function), indicating that confounding factors related to vascular complications were excluded.Only HDL-C, homocysteine concentrations, and white blood cell counts differed significantly between all four groups in both the discovery and validation cohorts.
All T2DM patients received combination therapy comprising a biguanide derivative (metformin) and sodiumglucose co-transporter type 2 inhibitors (SGLT-2is) or glucagon-like peptide-1 receptor agonists (GLP-1RAs).Participants diagnosed with CAD, either with or without T2DM, were treated with low-dose aspirin and did not receive any other antiplatelet medications.Additionally, approximately 90% of individuals with T2DM and all patients with CAD in the discovery and validation cohorts received statins, resulting in normal lipid parameters or, in some cases, only mild dyslipidemia.The frequency of hypertension was higher among patients with T2DM-CAD (100%) and CAD (100%) compared to those with T2DM (91.67% vs. 96.67%) in the discovery and validation phase of the study (p = 0.423 and p = 0.460).No significant intergroup differences were observed with respect to the antihypertensive treatment regimen used.
The demographic and clinical characteristics of patients enrolled in the discovery and validation cohorts are summarized in Table 2.

Functional annotation and pathway enrichment analysis
As a particular miRNA may act on several targets, the DIANA-miRPath v3.0 tool was used to identify the top 10 most significant KEGG pathways and GO annotations corresponding to gene targets for each of the 12 DE-miRNAs between patients with T2DM-CAD and T2DM.Subsequently, possible relationships between the selected up-and down-regulated miRNAs were visualized using Cytoscape v3.10.0 software.The interaction network functionally linked the 11 miRNAs together, leaving only one down-regulated hsa-miR-4706 without functional linkage (Fig. 2).
To further recognize the specific biological functions of DE-miRNA-regulated target genes, functional enrichment analysis was performed separately for upand down-regulated miRNA sets using the miRNet 2.0     3A).In the group of down-regulated miRNAs, the most prominent pathways were neurotrophin signaling pathway along with the signaling pathways of Wnt and mitogen-activated protein kinases (MAPK) and others associated with signal transduction via NGF, Rho GTPases and transforming growth factor-β (TGF-β) receptor complex.What is more, they have been found to regulate processes related to transcription and translation, the endocrine system, especially the insulin  signaling pathway, cellular community, including focal adhesion, and zinc ion binding (Fig. 3B).

Protein-protein interaction network construction
The PPI networks of the validated target genes of 12 DE-miRNAs were constructed separately for up-and downregulated miRNAs between the T2DM-CAD and T2DM groups (Additional file 1: Table S5).Using the STRING database in Cytoscape v3.10.0 software, we identified 319 nodes and 484 edges, with a PPI enrichment p-value of 5.3 × 10 −4 for up-regulated miRNAs, while a total of 1,710 nodes and 12,287 edges, with a PPI enrichment p-value < 1.0 × 10 −16 , were recognized for down-regulated miRNAs.
Next, a clustering analysis of all nodes and edges was performed using the MCODE algorithm in Cytoscape v3.10.0, and the 6 highest ranked clusters for both up-(clusters I-VI) and down-regulated (clusters VII-XII) miRNAs are shown in Fig. 4.
In addition, a functional enrichment analysis was conducted to gain insights into the biomolecular significance of each cluster's genes (Additional File 1: Table S6 and Table S7).Cluster I was associated with RNA metabolism, while clusters V, IX and XI were found to be functionally linked to multiple steps in transcription and translation processes.Clusters II and VII were predominantly involved in the cell cycle.As expected, DE-miRNA target genes, grouped into sub-networks III-IV, VI, VIII-XII, were strictly enriched in various pathways underlying the pathogenesis of diabetic atherosclerosis.Cluster III was connected with the Hippo-signaling pathway, cell-extracellular matrix-interactions and, together with cluster IX, was associated with focal adhesion.Cluster IV focused on the function of the phagophore assembly site membrane, which is the main structure involved in the induction of autophagy.Cluster VI members dealt with cation, mainly calcium channel complex and the regulation of cation transmembrane transport, whereas cluster X was functionally related to ion binding.Interestingly, cluster VIII was linked to the regulation of metabolic processes, vascular endothelial cell proliferation and migration, response to oxygen-containing compounds, leukocyte activation, platelet degranulation, and lipid homeostasis.The results of the functional enrichment analysis revealed that cluster XI was involved in signaling through the TGF-beta receptor complex by regulating the transcriptional activity of the SMAD2/SMAD3:SMAD4 heterotrimer, while cluster XII was associated with the cytoplasmic stress granule, which are dynamic cytoplasmic aggregates formed in response to cellular stress.

Hub gene identification
The PPI networks were then filtered using the cyto-Hubba plugin in Cytoscape v3.10.0 software to identify genes showing the high whole-network connectivity.For up-regulated miRNAs, the top 10 hub genes (Hub I) were Phosphatase and Tensin Homolog (PTEN), Protein Kinase C Alpha (PRKCA), Pyruvate Kinase M1/2 (PKM), CD4 Molecule (CD4), RNA Polymerase II, I and III Subunit E (POLR2E), Calreticulin (CALR), Fig. 3 Functional enrichment analysis for genes targeted by DE-miRNAs between T2DM-CAD and T2DM.The top 10 terms of the three Gene Ontology (GO) subcategories: Biological Process (GO: BP), Cellular Component (GO: CC), Molecular Function (GO: MF) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathways are presented separately for (A) up-regulated miRNAs and (B) down-regulated miRNAs.The terms are sorted for each category in descending order according to the adjusted p-values.p-value-EASE score for enrichment adjusted by the Benjamini-Hochberg correction for multiple hypothesis testing.Adjusted p-value < 0.05 was statistically significant.The number in brackets following the term names indicates the number of enriched target genes.The plot was generated using the ggplot2 3.3.0package in the R environment.DE-miRNA, differentially expressed miRNA; T2DM-CAD, type 2 diabetes mellitus with coronary artery disease; T2DM, type 2 diabetes mellitus Mitogen-Activated Protein Kinase Kinase 7 (MAP2K7), Sterol Regulatory Element Binding Transcription Factor 1 (SREBF1), U2 Small Nuclear RNA Auxiliary Factor 2 (U2AF2), and Solute Carrier Family 2 Member 1 (SLC2A1), as per the decreasing order of degree values (Fig. 5A).
Among the hub gene networks for up-and down-regulated miRNAs, PTEN and ACTB were designated as key hub genes with the highest node degrees of 72 and 119, respectively.Through the DE-miRNA-hub gene network construction, most of the hub genes were found to be potentially modulated by up-regulated hsa-miR-4505 and hsa-miR-4743-5p and down-regulated hsa-miR-3613-3p and hsa-miR-4668-5p (Table 4).
As listed in Table 5, Hub I genes were significantly enriched in response to stress, insulin, and oxygencontaining compound, as well as regulation of ErbB and hypoxia-inducible factor 1 (HIF-1) signaling pathways.Among the top 10 biological functions of Hub II genes, TGF-β and Hippo signaling pathways, and platelet activation were most associated with the development of diabetic atherosclerosis.Moreover, based on KEGG pathway analysis, several other T2DM-CAD-related processes, including leukocyte transendothelial migration, fluid shear stress and atherosclerosis, focal adhesion, along with apelin, Wnt, and chemokine signaling pathways, were found to enrich Hub II genes (Table 5).

Validation of differentially expressed miRNAs by RT-qPCR
The differential expression of five miRNAs (hsa-miR-3613-3p, hsa-miR-4505, hsa-miR-4668-5p, hsa-miR-4743-5p, hsa-miR-4750-3p) selected from miRNA profiling using microarrays based on the highest FC, AUC value and their potential biological significance in the pathogenesis of CAD in T2DM, was validated by RT-qPCR in a separate larger replication cohort comprising of 94 patients (T2DM with CAD group, n = 30; T2DM group, n = 30; CAD group, n = 16; control group, n = 18).The expression levels of miRNAs were normalized to the expression of U6 snRNA as an endogenous control.hsa-miR-4505, hsa-miR-4743-5p, and hsa-miR-4750-3p were constantly detected in all samples, but hsa-miR-3613-3p and hsa-miR-4668-5p could only be detected at very low levels by RT-qPCR, making their valid quantification impossible.Hence, these two miRNAs were excluded from subsequent analysis.
Similar to the discovery phase of the study, hsa-miR-4505 and hsa-miR-4743-5p were found to be significantly up-regulated, whereas hsa-miR-4750-3p was down-regulated in both T2DM-CAD group compared to the T2DM group and the T2DM-CAD group compared to the control group (all p < 0.001) (Fig. 6).The most increased miRNA expression was detected in hsa-miR-4505 with a 6.64-fold change, and a 3.50-fold decrease in expression level was found in hsa-miR-4750-3p in T2DM-CAD patients compared to T2DM individuals.

Evaluation of the diagnostic value of tested miRNAs in T2DM-CAD
ROC curve analysis was conducted to assess the diagnostic value of the three identified miRNAs as candidate biomarkers for CAD in T2DM.All three miRNAs, hsa-miR-4505, hsa-miR-4743-5p, and hsa-miR-4750-3p,  showed AUCs reaching 0.876, 0.860, and 0.833, respectively, indicating that they possess good diagnostic potential in discriminating between T2DM-CAD and T2DM patients (Fig. 7).In addition, we calculated the Youden Index to determine cut-off values.The optimal cut-off points (log 10 RQ) for hsa-miR-4505, hsa-miR-4743-5p, and hsa-miR-4750-3p were 0.22, 0.43, and − 0.54, respectively.The corresponding sensitivities and specificities to these points were evaluated and are presented in Additional file 1: Table S8.

Identification of the best combination of miRNAs for CAD detection in T2DM
To investigate the combination of miRNAs with the strongest classification power for discrimination between T2DM-CAD and T2DM subjects, logistic regression models were developed.Based on the results of the miRNA expression levels obtained during the validation phase of the study, four models were derived: model 1 was built on the expression of hsa-miR-4505 and hsa-miR-4743-5p as independent variables, model 2 included hsa-miR-4505 and hsa-miR-4750-3p, model 3 was based on the expression of hsa-miR-4743-5p and hsa-miR-4750-3p, while model 4 included all three variables.To assess the diagnostic potential of the constructed models, ROC curves were plotted (Fig. 8), and the basic parameters and common quality measures of the models are summarized in Additional File 1: Table S9.All miRNAbased signature revealed higher diagnostic values compared to the those for miRNAs used separately.The data indicated that the best classification accuracy for detecting CAD in T2DM was achieved by the three-miRNA

Discussion
T2DM and CAD constitute two closely inter-related clinical entities whose concurrent prevalence is steadily increasing [37].It is now well established that T2DM accelerates the onset of CAD, however, the pathogenetic molecular mechanisms underlying diabetic atherosclerosis are not fully understood [3].Despite substantial improvement in the diagnostic strategies in cardiovascular diseases, effective non-invasive screening methods for early detection of CAD in T2DM are still urgently needed [1,7,9].
In the present study, we investigated the potential role of circulating miRNAs as novel biomarkers for CAD in patients with T2DM.To find a unique miRNA profile in plasma samples from T2DM-CAD patients, we performed microarray analysis using a broad panel of 2,578 human miRNAs that have not been previously assessed in coexisting CAD and T2DM.Our results revealed that the expression of twelve miRNAs was significantly altered in T2DM-CAD patients compared to T2DM subjects.Moreover, upon the series of bioinformatic analyses, we uncovered the potential biological significance of DE-miRNA target genes in the pathogenesis of diabetic atherosclerosis.To the best of our knowledge, we have demonstrated for the first time a specific three-miRNA model based on the combined expression of hsa-miR-4505, hsa-miR-4743-5p, and hsa-miR-4750-3p that can effectively differentiate between T2DM patients with and without CAD.Based on the miRNA expression profiling with microarray platform, we detected a significant increase in the expression levels of three miRNAs (hsa-miR-4505, hsa-miR-4743-5p, hsa-miR-6846-5p), while nine miRNAs (hsa-miR-3613-3p, hsa-miR-4668-5p, hsa-miR-4706, hsa-miR-6511b-5p, hsa-miR-6750-5p, hsa-miR-4750-3p, hsa-miR-320e, hsa-miR-4717-3p, hsa-miR-7850-5p) were down-regulated in T2DM-CAD patients compared to T2DM individuals.Among them, hsa-miR-4505 and hsa-miR-3613-3p appeared to be the most dysregulated miRNAs in T2DM-CAD.So far, only limited studies have exploited the large-scale methods to provide comprehensive analysis of miRNA patterns in patients with T2DM-CAD [38][39][40].Leveraging the NanoString nCounter technology, Bielska et al. found that six miRNAs (miR-615-3p, miR-3147, miR-1224-5p, miR-5196-3p, miR-6732-3p, and miR-548b-3p) were overexpressed in the serum of T2DM patients with ischemic heart disease as compared to those with uncomplicated T2DM [38].In turn, Zhang et al. uncovered a set of 138 extracellular vesicle-carried DE-miRNAs in T2DM-CAD patients by using small RNA sequencing [39].Nevertheless, these miRNAs have only been evaluated for their utility in distinguishing between T2DM-CAD patients and healthy controls [39].Similar to prior reports, most of the identified miRNAs in our study, with the exception of hsa-miR-4505, hsa-miR-3613-3p and hsa-miR-320e, have not been previously investigated in CAD and/or T2DM and appear to be novel in the genomic control of T2DM complicated by CAD.
To date, it has been shown that miR-4505 was downregulated in patients with acute myocardial infarction, whereas miR-3613-3p and miR-320e were up-regulated in patients with stable CAD, but the pattern of changes was not consistent with those obtained in our study [41][42][43].These discrepancies in miRNA expression levels can be partially explained by the different sample types (plasma, serum, exosomes), the analytical techniques utilized for miRNA profiling, and differences in ethnicity and disease stages among the study populations.Nonetheless, similar to the current observation, a downregulation of the 5′-isomiR of miR-4717 was noticed in patients affected by acute myocardial infarction [44].Importantly, plasma miR-4717-5p was found to distinguish patients with ST-segment elevation and non-STsegment elevation myocardial infarction from controls with AUCs above 0.800, supporting its potential as a diagnostic biomarker for CAD [44].
Functional annotation analyses of our DE-miRNAs revealed an overrepresentation of genes in terms related to organelle, cellular nitrogen compound metabolic process, biosynthetic process, and ion binding.In addition, the targets of hsa-miR-4505, hsa-miR-4743-5p and hsa-miR-4750-3p were significantly enriched in pathways relevant to the initiation of atherogenesis, including fatty acid metabolism, the neurotrophin signaling pathway, and leukocyte transendothelial migration.
Metabolism of fatty acids engages several inter-connected processes, including fatty acid biosynthesis and degradation, which we found to be modulated by hsa-miR-4505, hsa-miR-4743-5p and hsa-miR-4750-3p [45].Disturbed fatty acid metabolism has been implicated in diabetes, as evidenced by elevated levels of circulating triglyceride-rich lipoproteins, saturated fatty acids, and increased fatty acid oxidation [46].Importantly, endothelial cells exposed to hyperglycemia and insulin resistance can rapidly switch from aerobic glycolysis as the predominant source of energy towards enhanced fatty acid β-oxidation [46,47].Since a significant amount of NADH is generated, the disturbances in redox homeostasis and increased production of harmful reactive oxygen species (ROS) occur [48,49].Moreover, excess free fatty acids lead to activation of protein kinase C (PKC), the nuclear factor-kappa B (NF-κB) pathway, and an increase in endothelial activation markers such as intercellular adhesion molecule-1 (ICAM-1) and vascular cell adhesion molecule-1 (VCAM-1) [50].Thus, the dysregulation of fatty acid metabolism triggers a myriad of detrimental effects, aggravating oxidative stress, endothelial permeability and inflammation, ultimately contributing to endothelial injury and atherogenesis [50].
In a recent study, and in agreement with our microarray data, up-regulation of miR-4505 was found to exacerbate endothelial injury, suggesting its potential as a novel proatherogenic molecular target [51].NF-κB signaling activated in response to pro-inflammatory stimuli in endothelial cells has been shown to induce miR-4505 expression and, by targeting heat shock protein A12B, enhance endothelial permeability by reducing cell-cell junctions, especially vascular endothelial (VE)-cadherin [51].Similarly, Liu et al. confirmed the atheroprone effect of miR-3613-3p down-regulation in heat-stressed endothelial cells, identifying the mitogen-activated protein kinase kinase kinase 2 (MAP3K2)/p38/caspase-3 pathway as its downstream target [52][53][54].As previously reported, we reinforced that hsa-miR-3613-3p can modulate the Hippo signaling pathway, consisting of a cascade of kinases, transcriptional coactivators and critical effectors such as Yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) [55].A variety of stress signals, including hypoxia, endoplasmic reticulum stress or heat stress, can trigger the Hippo signaling pathway [55].Interestingly, inactivation of Hippo/YAP signaling has been found to promote endothelial inflammation in a model of diabetes-accelerated atherosclerosis by enhancing monocyte-endothelial cell adhesion and production of pro-inflammatory cytokines [56].Additionally, Hu et al. stated that inhibition of the Hippo-YAP/TAZ/miR-496 pathway in oxidized-LDL (ox-LDL) treated endothelial cells leads to their dysfunction [57].Therefore, the exact role of the miR-3613-3p/ Hippo signaling pathway in the pathogenesis of diabetic atherosclerosis needs to be clarified.
A comprehensive bioinformatic analysis highlighted that pathways in cancers, calcium ion binding, cell-tocell communication, and actin cytoskeleton dynamics were among the most enriched terms for the set of upregulated miRNAs in our study.It is well recognized, that under hyperglycemic milieu, the inflamed endothelium changes its structural and functional properties, mainly through RhoA, a small Rho GTPase [58].Mechanistically, chemokine-induced activation of these intracellular signaling proteins increases cytoskeletal dynamics and contractility, and leukocyte polarization, enabling directional binding of endothelial adhesion molecules to their respective ligands on the leukocyte surfaces [59].Finally, calcium-dependent VE-cadherin-based adherens junctions are disrupted, disassembled, internalized and ultimately degraded, enhancing endothelial permeability and promoting leukocyte transendothelial migration, a key step in atherosclerotic plaque formation [58,60].Hence, we proposed hsa-miR-4505, hsa-miR-4743-5p and hsa-miR-6846-5p as novel potential atheromiRs, but in vivo and in vitro studies are needed to further investigate the role of these miRNAs in diabetic atherosclerosis.
Among the DE-miRNAs identified in T2DM-CAD, we discovered hsa-miR-4505, hsa-miR-4717-3p, and hsa-miR-7850-5p to be essential regulators of signaling by neurotrophins, predominantly NGF, through TRK receptors, which are highly expressed in endothelial cells and VSMCs [61].Emerging evidence indicates that the expression of NGF dramatically increases in response to vascular injury, providing the activation of several pathways, including MAPK/extracellular signal-regulated kinase (ERK), phosphatidylinositol 3-kinase (PI3K)-Akt, and the small GTPase RhoA, the latter of which promotes VSMC proliferation and migration [61].In contrast, Chaldakov et al. showed decreased levels of NGF in plasma samples from patients with metabolic syndrome and in the atherosclerotic coronary vessel wall, especially those with advanced lesions, confirming its involvement in the evolution of atherosclerotic plaques in T2DM [62,63].
PTEN acts as a stress-balancing molecule and is involved in the regulation of endoplasmic reticulum stress, endothelial cell apoptosis, inflammation, and VSMC proliferation and migration in response to ox-LDL-induced oxidative stress [64][65][66][67][68][69].Moreover, calcium-dependent and diacylglycerol-activated PKCα (PRKCA) under hyperglycemic settings has been found to promote ROS-producing enzymes and oxidative stress, the expression of endothelial adhesion molecules, proinflammatory cytokine production, and deteriorate nitric oxide-dependent vasodilation, indicating that it may be a mediator of the onset of diabetes-related atherosclerosis [70].Similarly, CALR, encoding calreticulin, an endoplasmic reticulum stress marker, and MAP kinase kinase 7 (MAP2K7), a key upstream transducer of stress-activated protein kinase, have been evidenced as vital contributors to endothelial dysfunction in an oxidative stress environment [71,72].In line with previous reports, the results of our functional enrichment analysis of PTEN, PRKCA, CALR, MAP2K7 argue for their involvement mainly in the regulation of the oxidative stress response in diabetes.Additionally, recent studies have suggested that monocytes and macrophages from patients with atherosclerotic CAD overuse glucose and overexpress several glycolysis-related genes, including PKM and SLC2A1, encoding two PKM isoforms (PKM1 and PKM2) and glucose transporter 1 (GLUT-1), respectively [73,74].PKM2 and GLUT-1 have been demonstrated to promote the aberrant differentiation of macrophages towards a proinflammatory M1 phenotype, enhance ox-LDL uptake, thereby fostering foam cell formation and exacerbating inflammation [73][74][75][76][77]. SREBF1, another identified hub gene in our study, plays a central role in regulating lipid biosynthesis and uptake, and its elevated levels were found during the progression of foam cell formation in ox-LDL-treated macrophages [78].Hence, PKM, SLC2A1 along with SREBF1 bridge metabolic and inflammatory dysfunctions in the pathogenesis of diabetic atherosclerosis.Fu et al. revealed the levels of CD4, which encodes a glycoprotein expressed on immune cells such as CD4 + T cells, are increased in atherosclerotic plaques, particularly those with advanced lesions [79].It is of note that CD4 effectively distinguished T2DM-CAD patients from controls [79].
As we mentioned above, β-actin, encoded by ACTB, is a core cytoskeleton protein that, as an effector of the small GTPase RhoA, regulates cell structure, proliferation and migration, contributing to vascular endothelial dysfunction and remodeling [80].The elevated expression of ACTB and MYC was also detected in epicardial adipose tissue samples from patients with CAD, implying their paracrine modulatory role in the development of atherosclerosis [81].Interestingly, activation of MYC, a key downstream target gene of the Wnt signaling pathway, has been observed at an early stage of atherosclerotic plaque evolution [82].Similarly, up-regulation of RhoA in a hyperglycemic state aggravates endothelial dysfunction and increases endothelial permeability through stress fiber formation, focal adhesion and cell contraction [83,84].Conversely, inhibition of RhoA, which is a downstream effector of Hippo/YAP/TAZ activation, exerts atheroprotective and anti-inflammatory effects [85].MYC, RHOA, SMAD2, known as TGF-β transducers, along with CDK1, have been identified as stimulators of VSMC proliferation, migration, apoptosis or ferroptosis, thereby exacerbating diabetic atherosclerosis [84,[86][87][88][89]. Deng et al. and Yang et al. postulated that disturbed fluid shear stress activates TGF-β receptors, amplifies SMAD2/SMAD3 activation, and thus the expression of pro-inflammatory genes [90,91].These observations may partially explain the increased expression of SMAD2 in endothelial cells derived from patients with CAD [92].In accordance with previous studies, our analysis of KEGG pathways confirmed the involvement of GNAQ in platelet activation [93].Notably, all of the identified GO terms, KEGG and Reactome pathways for hub genes of downregulated miRNAs are strongly related to the development and progression of diabetic atherosclerosis.
To verify the robustness of the microarray analysis, we further validated the expression of five selected DE-miRNAs using RT-qPCR in a larger independent cohort, confirming the overexpression of hsa-miR-4505, hsa-miR-4743-5p and the down-regulation of hsa-miR-4750-3p in T2DM-CAD patients with respect to T2DM individuals.Nevertheless, the expression of the highly down-regulated hsa-miR-3613-3p and hsa-miR-4668-5p could only be detected at very low levels, making their valid quantification impossible.This lack of concurrence between these two methods may be partially elucidated by the observed lower correlation for down-regulated genes in prior investigations due to the greater variability associated with decreased reaction efficiencies found at later cycles, where genes with low expression levels respond [94].
Evaluating the diagnostic accuracy of hsa-miR-4505, hsa-miR-4743-5p and hsa-miR-4750-3p, we showed that each of these miRNAs could be a good stand-alone discriminator of CAD in T2DM patients.These results are consistent with prior reports that have demonstrated comparable performance values for different single miRNAs in CAD associated with T2DM [15-18, 95, 96].Innovatively, we displayed that the model based on combining the expression of hsa-miR-4505, hsa-miR-4743-5p, hsa-miR-4750-3p yielded the highest diagnostic accuracy with an excellent AUC of 0.959 (95% CI, 0.914-1.000;p < 0.0001) and 100.00% sensitivity and 86.67% specificity in comparison with each miRNA individually and constructed panels of two miRNAs in detecting CAD in T2DM patients.Likewise, other authors suggested that the two miRNA signatures consisting of miR-3147, miR-615-3p and miR-9, miR-370, respectively, revealed better diagnostic value than these miRNAs considered individually in identifying T2DM patients at high risk of CAD [38,97].Although Al-Muhtaresh et al. also achieved good discriminatory efficacy for the combined miR-1 and miR-133 signature (AUC 0.752), this was lower than that for miR-1 alone (AUC 0.802), possibly due to only a marginally statistically significant AUC value for miR-133 [98].These findings provide the groundwork for future investigations, whereby a panel of miRNAs appears to be more beneficial in predicting the development of CAD in patients with T2DM rather than single miRNAs.

Limitations of the study
The present study also has some limitations.Firstly, due to the observational nature of the study, the potential impact of residual confounding factors cannot be ruled out.Nevertheless, our investigation was based on the rigorous methodological design and comprehensive data collection methods.Secondly, the sample lacked ethnic diversity, which limits the generalizability of the findings.Thirdly, although functional and pathway enrichment analyzes were performed to determine the potential biological significance of DE-miRNAs, further studies in in vitro and in vivo experimental models are needed to conclusively define the regulatory role of miRNA target genes in the pathogenesis of CAD in T2DM.Fourthly, the lack of inter-laboratory standardization of miRNA detection and the costs associated with miRNA analysis may partially restrict the translation of miRNA-based biomarkers into clinical practice.Finally, the sample size is limited, thus the results should be considered preliminary and verified in large-scale, multi-center studies to confirm the applicability of the three-miRNA signature as a diagnostic biomarker for CAD in patients with T2DM in the general population.

Conclusions
In summary, the present study revealed a unique profile of circulating plasma-derived miRNAs in T2DM patients with comorbid CAD using microarray technology.Integrated bioinformatics analysis showed that DE-miRNA target genes were mainly enriched in pathways related to regulation of cellular processes, gene expression, fatty acid metabolism, oxidative stress response, neurotrophin signaling pathway, and leukocyte transendothelial migration, known to be involved in the initiation of diabetic atherosclerosis.We independently confirmed that hsa-miR-4505 and hsa-miR-4743-5p are significantly overexpressed, while hsa-miR-4750-3p is down-regulated in T2DM-CAD patients compared to T2DM subjects and may serve as novel non-invasive molecular markers of diabetes-accelerated CAD.Innovatively, we developed a three-miRNA model based on the combined expression of hsa-miR-4505, hsa-miR-4743-5p, hsa-miR-4750-3p, which could effectively detect CAD in T2DM in the Caucasian population, better than each miRNA separately.
tool.The top 10 KEGG and Reactome pathways and GO terms categorized as the biological process (BP), cellular component (CC), molecular function (MF), potentially relevant to the development of CAD in T2DM, are presented in Fig.3.The gene targets of the up-regulated miRNAs were significantly overrepresented in terms related to pathways in cancers, signal transduction, mainly calcium, epidermal growth factor receptor (EGFR) and nerve growth factor (NGF) signaling pathways, metabolism, calcium ion and protein binding, cell-cell communication, cell adhesion, and cell junction organization (Fig.

Fig. 2
Fig. 2 Functional network of DE-miRNAs in T2DM-CAD along with their top 10 pathways and biological features.Down-regulated miRNAs are marked in yellow, while up-regulated miRNAs are marked in green.Analysis was performed using the DIANA-miRPath v3.0 tool.DE-miRNA, differentially expressed miRNA; T2DM-CAD, type 2 diabetes mellitus with coronary artery disease

Fig. 4
Fig. 4 Analysis of the PPI interaction network for DE-miRNA target genes in T2DM-CAD.The 6 top-scored clusters selected from the PPI network for both (A) up-regulated and (B) down-regulated miRNAs are presented along with the number of nodes, edges, and key gene (highlighted by enlarging) in each cluster.Nodes having different colors indicate different proteins encoded by the co-expressed genes.Gray lines connect proteins within the PPI sub-networks with darker colors and thicker lines indicating higher core PPI values.The analysis was performed using the STRING database with the MCODE plugin of Cytoscape v3.10.0.PPI, protein-protein interaction; DE-miRNA, differentially expressed miRNA; T2DM-CAD, type 2 diabetes mellitus with coronary artery disease

Fig. 5
Fig.5 The top 10 hub genes in PPI networks of DE-miRNA target genes.(A) Hub I for up-regulated miRNAs.(B) Hub II for down-regulated miRNAs.The node color changes gradually from red to yellow in descending order according to the BottleNeck ranking method.PTEN and ACTB were recognized as key hub genes for up-and down-regulated miRNAs, respectively.The analysis was performed using the cytoHubba plugin of Cytoscape v3.10.0.PPI, protein-protein interaction; DE-miRNA, differentially expressed miRNA

Table 5
Functional enrichment analysis of hub gene modules for up-regulated and down-regulated miRNAs in T2DM-CAD T2DM-CAD; type 2 diabetes mellitus with coronary artery disease; GO, Gene Ontology; BP, biological process; MF, molecular function; KEGG, Kyoto Encyclopedia of Genes and Genomes p-value and adjusted p-value < 0.05 were statistically significant